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Abstract 

We explore how heterogeneity in the intensity of interactions between people affects epidemic 
spreading. For that, we study the susceptible-infected-susceptible model on a complex network, 
where a link connecting individuals i and j is endowed with an infection rate (5ij = Xuiij pro- 
portional to the intensity of their contact u>ij, with a distribution P(wij) taken from face-to-face 
experiments analyzed in Cattuto et al. (PLoS ONE 5, ell596, 2010). We find an extremely 
slow decay of the fraction of infected individuals, for a wide range of the control parameter A. 
Using a distribution of width a we identify two large regions in the a — A space with anoma- 
lous behaviors, which are reminiscent of rare region effects (Griffiths phases) found in models 
with quenched disorder. We show that the slow approach to extinction is caused by isolated 
small groups of highly interacting individuals, which keep epidemic alive for very long times. 
A mean-field approximation and a percolation approach capture with very good accuracy the 
absorbing-active transition line for weak (small a) and strong (large a) disorder, respectively. 

PACS numbers: 89.75.Hc,05.40.-a,87.19.X-,02.50.Ey 
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Nowadays, much of the daily human activity is constantly being recorded by means 
of modern technologies such as Internet, GPS, mobile phones, blue-tooth and other elec- 
tronic devices. The gathering and analysis of large amount of activity data have allowed 
to explore statistical features on people's behavior. In particular, recent studies have 
revealed interesting properties about how humans interact with each other, either by hav- 
ing conversations l|, or sexual contacts Q, or by means of mobile devices Q or wireless 
communications |4j. Despite the fact that the type of interactions in these studies are 
quite different, they all found a common interaction pattern, that is, human contacts are 
very heterogeneous. This is consistent with a broad distribution of different magnitudes 
that quantify the timing of contacts, such as their duration, frequency and gaps. This 
diversity could eventually have an impact on propagation processes that involve human 
contact, such as the spreading of rumors or diseases. 

In this letter we explore epidemic spreading on a population with heterogeneous in- 
teraction intensities. We use a distribution of intensities extracted from the pattern of 
contacts between participants of a conference obtained in recent face-to-face (F2F) 
experiments Jl|. For the spreading process we use the Susceptible-Infected-Susceptible 
(SIS) dynamics on Erdos-Renyi (ER) networks with infection rates across links 
that are proportional to the intensity of encounters. 

As the infection rate increases, the SIS model 6| exhibits a transition, from an absorb- 
ing (disease-free) phase where the infection dies exponentially fast to an active (endemic) 
phase where the infection spreads over a large fraction of the population and becomes 
persistent. We find that the heterogeneity in the intensity of contacts introduces an 
intermediate absorbing region, in which the epidemic dies very slowly, as a stretched ex- 
ponential or a power law in time. We experiment with other rate distributions and show 
that this slow approach to epidemic extinction is caused by the presence of small clusters 
composed by links with high interaction time, which remain infected for very long times. 
We also discuss analogies with the effects observed in models with quenched disorder 8|. 

While our results are mainly concerned with the decay of the infection in the epidemic- 
free phase, some related models |9-11| have focused, instead, on the disease prevalence 
within the endemic phase. Other studies have introduced heterogeneity at the individual 



evel, by assigning power law intertime events 
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12L Il3| or node-dependent infection rates 



14|. In our model heterogeneity is at the interaction level, by means of link-dependent 
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infection rates. 

SIS dynamics with F2F disorder. In the SIS model [6j, each individual of a population 
can be either susceptible (healthy) or infected. Infected individuals transmit the disease to 
its susceptible neighbors in the network at a rate v and return to the susceptible state at 
a rate 7. The dynamics is controlled by the rescaled infection rate A = f/7. For A above 
a critical value A c , even a small initial fraction of infected nodes is able to propagate the 
disease through the entire network (active phase), while for A < A c the disease quickly dies 
out (absorbing phase), following an exponential decay in the number of infected nodes. 

This model describes disease spreading in an ideal population where transmission rates 
between individuals are all the same. However, in real populations we expect interactions 
to be heterogeneous, having a broad range of intensities, as recently measured by analyzing 
mobiles phone data |s| and by means of person-to-person experiments In order to 
explore how the behavior of the SIS model is affected by the heterogeneity of interactions, 
we run simulations of the dynamics on ER networks with infection rates distributed 
according to the weight distribution P[w) of F2F experiments (see Fig. [T]). In these 

experiments, participants of a three-day congress were asked to wear a Radio Frequency 
Identification device on their chest, so that when two persons were close and facing each 
other a relation of face-to-face proximity was registered. The weights w of Fig. Q] are 
defined as the total number of packets exchanged (or total contact time) between pairs of 
participants during the three days. 

We are assuming that infection rates are proportional to the total time individuals are 
in contact with each other, as the likelihood of transmission increases with exposure time 
-longer contacts imply a higher risk of infection-. Therefore, we assign an effective rate 
of infection = Xwij between two individuals i and j that are connected by a link of 
weight Wij, where A is a free parameter that acts as a transformation scale of contact 
intensities into infection rates. 

In Fig. [2] we show simulation results of the time evolution of the average density of 
infected individuals p, over many realizations of the SIS dynamics, starting from a config- 
uration where a small fraction of nodes have been randomly infected, and with infection 
rates following the distribution of Fig. [TJ All simulations in this article correspond, unless 
noted, to ER networks of mean degree (k) = 4 and N = 10 5 nodes. We found that, be- 
sides the typical behavior observed in the active and exponential phases of the clean SIS 
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FIG. 1: Probability distribution of face-to-face contacts intensities (weights w) of the 25th Chaos 
Communication Congress in Berlin, on a log-log scale. Intensity is defined as the total number 
of packages exchanged between two attendees, which is proportional to the contact duration. 
The large dispersion of the data reflects the large heterogeneity in the duration of contacts. 
Weights are rescaled to the interval [0, 1] for a better comparison with the theoretical distribution 
P(w) = 1/aw in the interval [e _a ,l], as shown in the inset for a = 6 (squares) and a = 3 
(diamonds). 

model (all infection rates are the same), there is an intermediate region between A = 0.25 
and 3.7 with very slow relaxation to the absorbing state. The region 2.0 < A < 3.7 is 
characterized by a power-law decay with a continuously varying exponent, while in the 
region 0.25 < A < 2.0 the decay is faster than a power law but slower than exponential 
(see the inset of Fig. |2]), and can be fitted by a stretched exponential. 

SIS with variable disorder strength. In order to understand this phenomenon we explore 
the dynamics for different distributions of weights. We assign to each link ij a weight 
Wij = e~ aTtj , where is a random number taken from a uniform distribution in the 
interval [0, 1], and a is a parameter that sets the range of in [e" a , 1]. This method 
generates a power-law distribution P(w) = 1/aw. The parameter a controls the width of 
the distribution, and measures the heterogeneity or strength of disorder. In the inset of 
Fig. [T]we plot P(w) for a = 6 and a = 3, which is intended to mimic the broad distribution 
of F2F contacts, even though the decay exponents are somewhat different. 

The distribution of infection rates is given by 

P(/3) = ^, with/Hn [Ae-<\A]. (1) 
Notice that when a — > we recover the clean model where = A for all ij. This kind 
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FIG. 2: Average density of infected nodes p vs time t, on a log-log scale, under the SIS dynam- 
ics on ER networks with (k) = 4 and TV" = 10 5 nodes. The infection rate distribution P(/3) 
corresponds to the weight distribution P(w) of Fig. [TJ with (3 = Xw, for values of the control 
parameter A = 3.9,3.7,3.6,3.5,3.4,3.3,3.0,2.5 (main plot) and A = 1.5,1.0,0.5,0.2 (inset), from 
top to bottom, p decays as a power-law for 2.0 < A < 3.7, as a stretched exponential for 
0.25 < A < 2.0, and as an exponential for A < 0.25, as shown in the inset on a linear-log scale. 



of disorder was already used in several works on complex networks [15|, [16[ , where it was 
bund that in the strong disorder limit a — > oo jl7|, magnitudes such as the current flow 



18| are highly affected, since links with very small weights are rarely traversed by the flow. 
In the same way, we expect that high-weight links facilitate the spreading of infections in 
our model, while low-weight links hinder the spreading. 

Numerical simulations show that the behavior of p under the theoretical disorder given 
by Eq. ([I]) (not shown) is very similar to the one observed in Fig. |2] for the F2F disorder. 
In the a — A phase diagram of Fig. [3] we summarize the different types of behaviors. 
Above the numerical transition line A™ (a) denoted by the red circles we find the active 
phase (white), where p reaches a stationary value larger than zero, and below we find the 
absorbing phase where p decays to zero. For a = we recover the clean transition point 
A° = A c (a = 0) = l/(k) = 0.25 of the clean model. The absorbing phase is divided into 
three regions. The exponential region (green), which appears for A < A°, characterized 



by the decay p 



~ e 



-at 



of the clean model, the weak effects region (yellow) where we 



observe an stretched exponential behavior p 



-at" 



(b < 1.0), and the strong effects 



region (orange), with a power law decay p ~ t~ 7 . Exponents a, b and 7 vary continuously 
with A and a. We note that the pure stretched exponential decay is strictly observed just 
above the line A = A°, while a smooth crossover to the pure power law behavior is found 
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FIG. 3: Phase diagram of the SIS model with infection rate distribution P(/3) = l/a/3, and (3 in 
[Ae _a ,A]. Colored regions correspond to the absorbing phase: orange and yellow for the strong 
and weak effects regions, and green for the exponential decay region. The dashed and solid lines 
are the MF [Eq. (|3j) ] and percolation [Eq. ([6])] approximations, respectively, for the transition 
between the active and absorbing phases. 

in the intermediate region. 

In summary, the heterogeneity in infection rates induces a large new region inside the 
absorbing phase, in which the temporal evolution of p exhibits an anomalous slow decay. 
This is caused by the presence of exponentially small isolated regions in the network 
where the system is locally active, that is, with infection rates > A°, which are able to 
sustain the activity for very long times. To check this, we calculated the size distribution 
of clusters composed only by infected nodes nj(s), at a fixed large time. Results are shown 
in Fig. |H Inside the weak and strong effects regions (A = 0.75 and 1.4), nj(s) is close to 
an exponential, and the size of the largest cluster s max is much smaller than the network 
size N = 10 5 (see inset of Fig. 0J). Also, the values 0.27 and 0.70 of the average infection 
rates inside these clusters for A = 0.75 and 1.4, respectively, show that the long-time 
activity is located inside active clusters, in which the average rate of infection (/?) > A°. 
For comparison, in the active phase (A = 2.0) is 7600 < s max < 9500, indicating the 
spreading of the disease over a large fraction of the network. 

Similar anomalous behaviors are found in models with disorder, given rise to the so- 
caUed Griffiths phases (GP) ^MR. The cognation of exponent rar e reg ions 
in space that survive for exponentially long times results in an overall slowing down of 
the dynamics, as we show below. The long-time contribution of active clusters to p is 
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FIG. 4: Main: cluster size distribution of infected nodes ni for a = 6 and values of A inside three 
different regions of Fig. [3j A = 0.75 (circles) and 1.4 (squares) in the weak and strong effects 
regions, and 2.0 (diamonds) in the active phase. Distributions correspond to snapshots of the 
network at fixed large times. The average infection rates inside each cluster are (/3) ~ 0.27, 0.70 
and 0.75 for A = 0.75, 1.4 and 2.0, respectively. Insets: size distribution of the largest cluster 
n max showing the appearance of a large component in the active phase. 



estimated as 



P 



where P(s) 



~ e 



-ps 



dssP(s) e~ tMs \ 



(2) 



2l| is the fraction of active clusters of size s and r (s) is the mean decay 
time of those clusters. By doing a saddle-point analysis, and using the finite-size scaling 
r(s) ~ e cs , one arrives to the power-law decay p ~ £~p/ c (with logarithmic corrections) 
observed in the strong effects region of Fig. [21 The size of active clusters is of order two 
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for A just above A°, leading to exponentially weak effects of the form p ~ 

Mean-field and percolation approaches. Within a mean-field (MF) approximation, p 
evolves according to p — —p + X{w){k)p(l — p), where \{w) = A(l —e~ a )/a is the average 
infection rate, and (k)(l — p) is the average number of susceptible neighbors of an infected 
node. The stationary solutions p = and p = (A — A^' IF ) /A correspond to the absorbing 
and active phases, respectively, with the transition point at 



AHa) 



(3) 



(fc)(l -e~ a Y 

Impurities, in the form of low-weight links, locally reduce infection rates, thus the transi- 
tion happens at a value Af F (a) > A°. 

Expression (j^J) (dashed line in Fig. [3]) is a very good estimate of A° um (a) for a < 14 (weak 
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disorder), but systematic deviations appear as a increases. Discrepancies arise because 
MF assumes that all links can spread the disease but, when a is large, a fraction of links 
have such small rates (inactive links) that infection never passes through them during the 
epidemic's life time, and thus the effective network for the spreading dynamics is diluted 
respect to the original network. When dilution is large enough the effective network gets 
fragmented into many small disconnected components and, as the disease cannot spread 
out of these components the active state is never reached. Therefore, the active-absorbing 
transition point for a large (strong disorder) corresponds to the percolation threshold. This 
occurs when the fraction of inactive nodes q (nodes attached only to inactive links) exceeds 
the critical value q c . For ER networks with Poissonian degree distribution = e~^^r 
is 

q = J2 l i p k = e~ {k)(1 ~ h \ where (4) 

k 

is the fraction of inactive links, and /3 m is the largest infection rate that does not allow to 
transmit the disease. At the percolation threshold (k) = (1 — f^) -1 in the iV — > oo limit, 
thus q c = exp [(1 — q,f)~ l In (/3 m e a / A£ orc ) /a — 1], from where the percolation transition line 
is 

= Pm • (6) 

Using q c = 0.7443 for a network of size N = 10 5 22], expression fl6]) with f3 n 



is in excellent agreement with A° um (a) for a > 14 (solid line in Fig. [3]). The value of 
(3 m is estimated from the crossover conditions A c MF (a) = A perc (a) and d\ c MF (a)/da = 
<9A c perc (a)/<9a between the MF and percolation lines at the weak-strong disorder crossing 
point a = a*. We obtain (3 m = — [elngc]" 1 and a* = —[(1 — q c ) lngj" 1 with a* ~ 13.243 
for the network used here. 

Summary and Conclusions. In summary, the heterogeneity in the intensity of contacts 
between individuals induces a regime with extremely slow (power-law or stretched expo- 
nential) relaxation to epidemic extinction, akin to the slowing down found in systems with 
quenched disorder. This effect is very robust, as it was observed using an empirical distri- 
bution of contact durations in F2F experiments, as well as a theoretical distribution with 
variable width. While temporal heterogeneity, causality and bursty activity was found 



to hinder spreading [ill, [13], we showed here that spatial heterogeneity has the counter- 
balanced effect, making epidemic more persistence by slowing down its extinction. Once 
a group of highly interacting individuals gets infected, is able to continuously reinfect 
each other at a high rate, keeping the infection inside the group for very long times. Our 
founding can be used to design efficient mitigation strategies for the disease. For instance, 
moderating the activity of highly interacting people could dramatically speed up the final 
stage of the epidemic. 

This work was financially supported by UNMdP and FONCyT (Pict 0293/2008). The 
authors thank Lucas D. Valdez for useful comments and discussions. 
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